Conservation of Angular Momentum in a Flux Qubit 
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Oscillations of superconducting current between clockwise and counterclockwise directions in a 
flux qubit do not conserve the angular momentum of the qubit. To compensate for this effect the 
solid containing the qubit must oscillate in unison with the current. This requires entanglement of 
quantum states of the qubit with quantum states of a macroscopic body. The question then arises 
whether slow decoherence of quantum oscillations of the current is consistent with fast decoherence 
of quantum states of a macroscopic solid. This problem is analyzed within an exactly solvable 
quantum model of a qubit embedded in an absolutely rigid solid and for the elastic model that 
conserves the total angular momentum. We show that while the quantum state of a flux qubit 
is, in general, a mixture of a large number of rotational states, slow decoherence is permitted if 
the system is macroscopically large. Practical implications of entanglement of qubit states with 
mechanical rotations are discussed. 

PACS numbers: 74.50.+r, 85.25.-j, 03.67.-a 



I. INTRODUCTION 



Flux qubits are formed by quantum superposition of 
current states in a superconducting loop interrupted by 
one or more Josephson junctions 1,2 . Quantum mechan- 
ics of such a qubit is described by a double-well poten- 
tial, similar to the textbook example^ of the ammonia 
molecule, NH3. In the latter example the tunneling be- 
tween states corresponding to the N-atom located to the 
left (L) or to the right (R) of the H 3 triangle creates 
quantum superposition of the \L) and \R) states, with 
the lowest energy doublet given by \L) ± \R), If one 
prepares the molecule in, e.g., the \L) state, the quan- 
tum mechanical average of the position of the N-atom 
oscillates harmonically between left and right at the fre- 
quency oj — A/h, where A is the energy splitting of the 
doublet. In the rigorous formulation of this problem the 
N-atom and the H3 triangle co-tunnel in such a manner 
that the position of the center of mass of the four atoms 
is preserved, thus conserving the linear momentum. 

In the simplest formulation of the flux qubit problem 
the role of left and right is played by clockwise and coun- 
terclockwise directions of the current. Typical values of 
the angular momentum associated with the current range 
from a few hundred h for a submicron SQUID loop^, to 
I0 5 h for a micron-size loop 5 , to 10 10 h for larger SQUIDs 6 . 
To conserve the angular momentum the tunneling of the 
current between clockwise and counterclockwise direc- 
tions must be accompanied by quantum transitions be- 
tween mechanical clockwise and counterclockwise rota- 
tions of the body containing the flux qubit. This creates 
a controversy^. Indeed, the co-tunneling of the super- 
conducting current and mechanical rotation needed to 
conserve the angular momentum requires entanglement 
of quantum states of the flux qubit with quantum states 
of a macroscopic body. In any reasonable experiment the 
phase of the wave function of the equipment containing 
the flux qubit must be destroyed instantaneously. Then 



how can the flux qubit preserve coherence on a measur- 
able time scale? This paper is devoted to the detailed 
analysis of the entanglement of current states with me- 
chanical rotations and its implications for superconduct- 
ing qubits. 

Within an exactly solvable model of a flux qubit em- 
bedded in an absolutely rigid rotator we obtain entangled 
eigenstates of the system and their dependence on the to- 
tal angular momentum J. When the system is prepared 
in the state with a certain direction of the superconduct- 
ing current, this state is, in general, a quantum mixture 
of many rotational states of the body. However, only 
tunnel splitttings Aj of the states belonging to the same 
J contribute to the oscillations of the superconducting 
current. We show that decoherence resulting from the 
broad statistical distribution over J is small as long as 
the body containing the qubit is macroscopically large. 
Thus, contrary to what one might think, the macroscop- 
icity of the body that is entangled with the qubit, is in 
fact required for low decoherence. We then study deco- 
herence of a flux qubit due to torques generated by the 
oscillating current in the elastic solid and show how de- 
coherence rates obtained within the two models match. 
Among other problems we discuss renormalization of the 
tunnel splitting by the elastic environment and superra- 
diant relaxation in a system of closely packed qubits. 

The paper is structured as follows. Exactly solvable 
quantum model of a flux qubit interacting with rotations 
of a rigid body is studied in [II] Quantum states of the 
qubit entangled with rotations of the body are obtained 
in Section III Al Section III Bl is devoted to decoherence 
due to rotational excitations of the body. Elastic en- 
vironment is considered in Section IIII1 The model that 
conserves the total angular momentum is formulated in 
Section IIII Al Section IIII Bl discusses decoherence of the 
flux qubit by internal torques. Renormalization of the 
tunnel splitting by the elastic environment is computed 
in Section llll CI Section llVl contains numerical estimates, 
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discussion of various effects originating from conservation 
of angular momentum, alternative interpretations of the 
results, and final conclusions. 



II. RIGID BODY 

A. Rotational states of a flux qubit 

First, we consider the tunnel-split states of a flux qubit 
and ignore conservation of the angular momentum. Let 
the lowest-energy doublet of a flux qubit be 



1 



V2 



(lt>±U», 



(1) 



where | f) and | \.) are the eigenstates of the operator of 
the angular momentum of the electronic current inside 
superconducting loop l z , 

lit) = l|t) 

L\i) = (2) 
Eigenfunctions ^± satisfy 

H^± = E±V± (3) 
with H being the Hamiltonian of the qubit and 

E- - E + = A (4) 

being the tunnel splitting. It is convenient to describe 
such a two-state system by a pseudospin 1/2. Compo- 
nents of the corresponding Pauli operator cr are 

^ - U)(tl + lt)UI 

<?y = iU)(t|-<|t)U| 

a z = it)(ti-u)ai. 

The projection of H onto | f) and | J.) states is 



(5) 



(m|i/|n) |m) (n| 



(6) 



According to Eq. (TTJ), 



*-) 



(7) 



It is now easy to see from Eq. (|3]) that 

(t t) = a \m i> = o 
ai#it> = <ti#u> = -A/2. 



(8) 



With the help of these relations one obtains from equa- 
tions © and © 



-{A/2)a x . 



(9) 



The general form of the wave function of our two-state 
system is 



2ft) 



(10) 



with | C | 2 + |C + | 2 = 1. If one imposes the initial condi- 
tion *(0) = | t), then 



and (l z ) = 1{<J Z ), with 



(a : )=«>(/)K|vT/(/))=<-«,s(^ 



(11) 



(12) 



This equation describes harmonic oscillations of the su- 
perconducting current at the frequency A/H between 
clockwise and counterclockwise directions. Another way 
to obtain this result is to use the equivalence^ of the 
Schrodinger equation for spin one-half to the precession 
equation for the expectations value of <x, 



dt\2/ \ Sa- 




which gives 



dt 
d I 



A, 



(14) 



The last two equations give Eq. (|Tj 

We shall account now for mechanical rotations of the 
body containing the flux qubit. In this Section we shall 
deal with an absolutely rigid body that can only rotate 
as a whole. As we shall see, this problem contains all of 
the components needed to understand the effects of en- 
tanglement required by the conservation of the angular 
momentum. Rotation by the angle <fi about the quanti- 
zation axis Z transforms the Hamiltonian of the qubit 
into 



(15) 



Noticing that the operator of the angular momentum of 
the superconducting current, l z (that is chosen in units of 
H), commutes with <f> it is easy to project this Hamiltonian 
onto | f) and | 1) . Simple calculation yields the following 
generalization of Eq. (|9]): 

H' a = Y, HH'\n)\m)(n\ 

m,n— ^,4- 

A r _OJ1At 



= "2 * 



„2ift 



[cos(2lcj))(T x + sin(2/0)cTj,] 



(16) 
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where a± — \{u x ± i&y)- 

To develop a rigorous formulation of the problem let us 
first assume that the body with the qubit is an isolated 
system in a pure quantum state described by a single 
wave function. The full Hamiltonian of the system is 

H = - j [a x cos(2^) + Cy sin(2ty)] , (17) 

where L z = —i(d/d(j)) and I = I z is the moment of inertia 
of the body for rotation about the quantization axis. It is 
easy to check that this Hamiltonian commutes with the 
operator of the total angular momentum, 



Jz — L z + l z 



■ d , 

~%~r~ + 1<J Z 



(18) 



Consequently, the eigenstates of (fTTf must be entangled 
states of l z and L z that are eigenstates of the total an- 
gular momentum J z : 



*.7±) = 



c 



J± 



\t) l ®\J-l) L ±^-\l) l ®\J+l) 1 



with J = J z . Simple calculation gives 



\ 



1 ± 



A 2 / 2 
A(M)' 2 (h.])' 2 



(19) 



(20) 



and 



E 



J± 



(hi) 2 (hj) 2 
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± 



A2 
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[hi) 2 {hjf 
i 2 



(21) 



for the energy levels. Here ± corresponds to =F in Eq. 
CEO) and J = 0,±l,±2,.... 

Alternatively, the same results can be obtained in the 
coordinate frame attached to the current loop. In this 
case one starts with the Hamiltonian 



fr {hL z f A 
11, = — (7 T 

21 2 



(hJ z - hi z f 

21 



Its eigenfunctions are 
1 



|*,7±>r 



V2 



(c J± |t)i±c JT U)0®|J) 



(22) 



(23) 



while eigenvalues are given by Eq. (|2"Tj) . The two coordi- 
nate frames are related by unitary transformation. 



B. Decoherence from rotations 

Any real macroscopic system should have some dis- 
tribution over J. According to Eq. (|2"Tj) . at large / the 
energies of the states corresponding to different J can be 
very close. Consequently, a macroscopically large num- 
ber of different J-states should contribute to the expec- 
tation value of any physical quantity. Since the phases 



of such states can differ significantly, the question then 
arises how the coherence of the flux qubit is influenced 
by this effect. Rigorous answer to this question is given 
below. 

To study decoherence, one should prepare the system 
in a state with a certain direction of l z , e.g. l z = +1, and 
study how (l z ) would depend on time. Naturally, the 
initial state should be obtained by subjecting the system 
to a strong bias field in the direction of the magnetic 
moment of the current loop. Adding the term —t\W<j z 
to the Hamiltonian, it is easy to work out the energy 
levels of the biased states: 



E 



J± 




W (M)(hJ) 
~2 + 1 



(24) 

For a large positive bias the states corresponding to the 
plus sign in the above equation have too high energies 
and can be ignored. In this limit the relevant energies, 
up to a constant, are 



Ej = E, 



h 2 (J-l) 2 [hL z f 



21 



21 



(25) 



in accordance with the expectation that they must be 
the energies of the rotational states of the body. To make 
sure that the system is magnetized in the direction of the 
field, that is l z = +/, it must be put in contact with a 
thermal bath at temperature T . This provides thermal 
distribution over Ej with probabilities given by 



p 1 f E J 



z = ^2 exp 



Ej 



(26) 

If at t = the field is removed and the system is isolated 
from the bath, it will be a mixture of J-states, 



\*Jl)o = \l)®\J~l), 



(27) 



with the probability of each J determined by Eq. (|26p . 
Time evolution of each J-state is provided by 



l*J*> 



Cj+ 
V2 



|*J+)e 



-iE J+ t/h 



+ 



V2 



Vj-)e 



-iE,,-t/h 



(28) 

Consequently, the time dependence of the expectation 
value of l z — lo z is determined by 



(29) 



Using the relations 



2 

J+ 



c 2 



J— J 



(* J+ |o-,|*j_) = Cj+Cj- (30) 
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one obtains 



ft 



ft 



-ft 



cos 



A.. 



(31) 



where 



E J+ -Ej-=AJl + p 



2(hi)(hJ) 
IA 



(32) 

Notice that only the energy splitting between states be- 
longing to the same J, separated by Aj, contribute 
to (<J Z ). For a given J / oscillations of the super- 
conducting current occur between (l z ) — I and (l z ) — 



Kft 



!)/(/? j + 1) as compared to the oscillations be- 



tween ±1 for J = ((3, j = 0). 

Formally, at T — 0, only the non-rotating state with 
J = I contributes to the sum in Eq. pip , providing 



ft 



1 



1 



ft 



l+ft 



A, 



(33) 



where Pi equals /3j at J = I. For a macroscopic body with 
a large moment of inertia Pi <C 1, so that the difference 
between Eq. (TT21 and Eq. (|3"3"]l is very small. The absence 
of decoherence at T — is related to the fact that the 
system is in a pure J-state. 

At T ^ rotations of a macroscopic body must be 
distributed over a macroscopically large number of J 3> /. 
Consequently, one can replace J — I in Eq. (|25|) with J 
and replace summation in Eqs. pip . (|26p by integration 
over J. This gives Z = \J2-KlksT jfi. Expectation value 
of o z depends on time through (A/h)t, 



1 



d^j exp 



r J- 



and is determined by a single parameter, 



(34) 



(35) 



Note that /3; = 2(M) 2 /(IA) contains a macroscopically 
large number I in the denominator. This provides 



Pi < /3t < 1 



(36) 



For a macroscopic body it is small due to the smallness 
of Pi . In this limit the time dependence of the oscillating 
term in Eq. (|34|) can be computed exactly: 



Re 



j(A/h)t 



y/l - i<7oo(A/ H)t 



(38) 



One can see that th e amplitude of quantum oscillations 
is decreasing as 1/ ' y/ '<Joo(A/h)t. Thus, the effective de- 
coherence rate due to the entanglement of the flux qubit 
with rotations of the rigid body is 



n n 



AM 2 fk B T 



(39) 



Notice that slow, decay of coherent oscillations 

given by Eq. (|38|) is a consequence of the absolute rigidity 
of the body. 

Proportionality of T r to 1/ 1 illustrates our point that, 
contrary to the naive picture that one might havei, the 
entanglement of a flux qubit with rotations of a macro- 
scopic body, dictated by the conservation of angular 
momentum, does not necessarily result in a strong de- 
coherence. This comes as a consequence of the selec- 
tion rule: According to Eq. (f3"Tj) only tunnel splittings, 
A j = Ej + — Ej-. of the states (|2~T1) belonging to the 
same J contribute to (cr z ). For a macroscopic body, all 
Aj are very close, thus providing low decoherence. 



III. ELASTIC BODY 
A. Flux qubit in the elastic environment 

Realistically, the body containing a flux qubit is not 
absolutely rigid. During half-period of oscillations of the 
superconducting current the elastic stress generated by 
the changing angular momentum of the current may only 
extend as far as half- wavelength, A/2 = irhvt/A, of the 
transverse sound of frequency A/h and speed Vt ■ We shall 
assume that this distance is greater than the size of the 
current loop. For, e.g., a micron-size loop this condition 
would be typically fulfilled for A/H < 10GHz. It allows 
one to treat the flux qubit as a point source of the elastic 
stress, without considering interactions of segments of 
the current loop with the elastic environment. 

Now the rotation angle (f> that appears in the previous 
section is determined to the elastic twist ^ 



for any reasonable values of I, A, and T. Since the 
main contribution to the integral in Eq. (J34J) comes from 
Pj ~ Pt 3> Pi, the overwhelming majority of J con- 
tributing to the integral satisfy J 3> i in accordance with 
our assumption. 

From Eq. (f3~4"|) the asymptotic value of (a z ) is 



= lim (a z 

t— yoo 



= t;Pt - 2/3, 



knT 



(37) 



[V x u ], , 



(40) 



where u is the phonon displacement field at the location 
of the flux qubit r = 0. Conventional quantization of 
phonons gives 



2pV 



E 

kA 



[ik x e kA ]. 



-kA , 



(41) 
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where a kA , a k A are operators of creation and annihilation 
of phonons of wave-vector k and polarization A, e kA are 
unit vectors of polarization, w k A = v t k is the phonon 
frequency, p is the mass density of the solid and V is its 
volume. Since we limit our consideration to elastic twists, 
only the two transverse polarizations of sound contribute 
to Eq. gl]). 

Elastic Hamiltonian that replaces Hamiltonian (fTTj) of 
the rigid-body approximation is 



<t + exp 



kA 



kA «kA — Ob 



<r_ exp 



kA 



kA I OkA — O k A 



where 



kA 



i h [k x e kA ], 

2/3^ ^kA 



(42) 



(43) 



Validity of this approximation relies on the fact that an- 
gular velocity of the local rotation, = d<fi/dt, is always 
small compared to the frequency of sound lo. Indeed, 
noticing that according to Eq. (|40]) ft ~ uiku we see that 
fl <C uj coincides with the condition of validity of the 
elastic theory: 

Unitary transformation H r = U~ 1 HU with 



U = exp 



kA I akA — a kA 



kA 



(44) 



transforms Hamiltonian (021 into 



L kA 



- A 
U--a x 



kA 



a kA«kA - -2~?kA ^kA + < A 



(45) 



where an insignificant constant has been omitted. In the 
transition from the first to the second line of Eq. (|45|) we 
have used properties of the displacement operator, 



L)- 1 (a)aD(a) =a + a, b~ 1 (a)a)b(a) = a 1 * + , 



with 



akA 



kA 



(46) 



(47) 



Eq. (|45p shows that from mathematical point of view 
the problem formulated in this Section is a variance of 
spin-boson problem 10 . While some important theorems 
have been proved for this problem in recent years (see, 



e.g., Ref . [Til and references therein), its exact eigenstates 
are unknown. This prevents us from developing rigorous 
mathematical approach to decoherence along the lines of 
the previous Section. From a physical point of view, the 
attractiveness of our variance of the spin-boson model is 
in the absence of free parameters. The boson field in our 
case is the phonon displacement field. Its coupling to the 
flux qubit (described by spin 1/2) is completely deter- 
mined by the conservation of total angular momentum. 
In what follows, we will use an approximation based upon 
observation that local twists of the elastic solid due to os- 
cillations of the superconducting current in a flux qubit 
must be very small. Within this approximation we will 
describe transverse phonons by a classical displacement 
field u(r, t), satisfying V • u = 0. 

Expanding Hamiltonian (|4"2")l to the lowest power on 
the elastic twist and replacing operators by their classical 
expectation values, one obtans 



H = H 



E~ 



A A 

■lf<7x--z-l(r y 



d 3 rS(r) 



du x 
dy 



dUy 

dx 



(48) 



where He is the Hamiltonian of free rotations, 

2 



1 



d 3 rpv 2 



du a dup Y 
drp dr a J 



(49) 



The dynamical equation for the displacement field is 
d 2 u a _ da a p 



dt 2 



(50) 



where a a p = 8H/5e a p is the stress tensor and e a p 
dua/drp is the strain tensor. This gives 



d 2 u x 
dt 2 

d 2 u v 



v 2 S7 2 u a 



A z d v ^ 
"2 v dy ^ 



A, d 
— In.. 
2 



K-J(r). 



(51) 



(52) 



The above equations should be solved together with the 
Landau-Lifshitz equation for cr: 



h da 



SH 



2 dt ° ' Scr 



which gives 



^ = -a z Al J d 3 r6(r) 



du x du y 
dy dx 



dt 



cr.A 



(53) 

(54) 
(55) 



d(7 f 

= ~<J y A + cr x Al / d 3 rS(r) 



dy dx 



It is easy to see that Eqs. (j5"4"l) . and (|5oT) preserve 
the length of cr: a 2 x + a 2 + a 2 = 1. 
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First, let us show that, in accordance with our general 
line of reasoning, the above equations conserve the Z- 
component of the total angular momentum, 



J? = hla z 



(57) 



Here L z is the ^-component of the mechanical angular 
momentum. Its time derivative equals the Z-component 
of the total mechanical torque, K z , acting on the body. In 
the absence of the external torque applied to the surface 
of the body, K z is given by^ 



K. 



= J d 3 r (a yx 



(58) 



Conventional elastic theory postulates no internal 
torques, in which case the stress tensor would be symmet- 
ric and K z would be zero. Situation changes when there 
are transitions between angular momentum states of a 
microscopic object inside the body, such as, e.g., a flux 
qubit. In this case the stress tensor is non-symmetric, 
yielding 



dL z 
dt 



(59) 



where a' af3 — 8H lrA /8e a p is the part of the stress tensor 
related to the interaction of the flux qubit with the elas- 
tic environment, Hi nt . The latter is given by the second 
term in Eq. (1421) . To prove conservation of the total an- 
gular momentum one needs to write this term with the 
accuracy to second-order terms on the elastic twists: 



This gives 



I . 



; . 8U X 8Uy 

4 rf(r ' 



(60) 



yx 



i i 2 r 

-Aa y S(r) + -Aa x S(r) / d 3 rS{r) 
Aa y S(r) - l -Aa x S(r) J d 3 rS(r) 



du x _ duy 
dy dx 
du x du v 
dy dx 



(61) 



so that 



dJz 
dt 



= hi 



da z 
~~dT 



l 2 Aa, 



dL 
~dt 

d 3 rS(r) 



m 



da z 

dt 
du x 

dy 



lAtJy 
8Uy 

dx 



(62) 



It is now easy to see that condition dJ z /dt = coincides 
with one of the equations of motion, Eq. <\5t 



B. Decoherence from internal torques 

At u = equations (I5"5"j) and (|56p would describe coher- 
ent precession of cr about the X-axis, with a x = const, 



a z oc cos(tA/fi), and a y cx svnitAjK). Conservation of 
angular momentum makes the flux qubit wiggle mechan- 
ically when the current oscillates between clockwise and 
counterclockwise. Consequently, it becomes a source of 
sound, as can be seen from Eqs. (IBT]) and (|52|) . Let us 
linearize all equations of motion around o~ x = 1, u = 0, 
with small <T y ,g(t) oc e~ %u)t and 



u Xj y(r,t) cx e" 



d 3 k 



r u x , y (k) . (63) 



Writing 8(r) as J e' k r one obtains from Eqs. ([ST|) 
and H2J 



ikyOy IA ik x a y 

^( k ) = -— ,. 2 ..2 — » u y( k ) = 7T: 



2p k 2 v\ — uj 2 



2p k 2 v 2 — uj 2 



( 64 ) 

where k 2 — k 2 +k 2 + k 2 . Substitution into Eqs. ([55]) and 
(l56l results in 



h 2 uj 2 = A 1 — 



l 2 A 



d 3 k 



2p J (2tt) 3 k 2 v 2 



(65) 



The integral in this equation should be computed in the 
complex plane with account of a small imaginary part of 



d 3 k 



(2tt) 3 k 2 v 2 - uj 2 
This gives 



1 

3^ 



k 4 dk 



k 2 v 2 — , 



2 2 W 2 



= A 1 - 



6npvf 



that is, 



where 



A T 

n 



r n = 



l 2 A 5 



12irh 4 pvf 



(66) 



(67) 



(68) 



(69) 



is the T = rate of the decay of the coherent precession of 
cr. This result is in full agreement with the decoherence 
rate computed with the help of the Fermi golden rule by 
considering spontaneous quantum transition from the ex- 
cited state — | — I)) to the ground state + | — I)) 
with the radiation of a phonon of energy Ai2. Its general- 
ization to finite temperature is T e = To coth[A/(2/csT)]. 
At k B T > A it gives r e oc T as in Eq. ([39]) obtained for 
the rigid body. Comparison of the decoherence provided 
by the two models will be done in Section IIVI 

As is clear from the derivation, the above result corre- 
sponds to the decoherence of a weakly excited state of the 
flux qubit. Our method, however, permits study of deco- 
herence of the state prepared with u = and arbitrary 
a z (including a z = 1) at t = 0. Dynamics of the vector cr 
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consists of fast precession about the X-axis and slow re- 
laxation towards the energy minimum that according to 



Eq. (|48j) corresponds to cr x = 1, a ViZ = 0. It is accompa- 
nied by radiation of sound due to the torque acting on the 
flux qubit from the oscillating current. Noticing that the 
space-time Fourier transform of the displacement gener- 
ated by the torque, u(k, w), and the time Fourier trans- 
form, cr(ui), of <r(t) are always related by Eqs. ([54]) due 
to the linearity of Eqs. (fBTj) and (|52j) . one can transform 
the integral in Eqs. (|54| and (|56|) as 



d 3 r( J (r) ,^ dUy 



I A f duj n t 



dy dx J 6irvf J 2tt 

(70) 

To the first approximation, fast-precessing and slowly- 
relaxing solution of Eqs. ([55]) and ([56]) that satisfies a 2 + 
°l + °l = 1 is 



a y (t) = y/l- (a x ) 2 sin ^t 

o z {t) = v/1- (^) 2 cos (~t\ , (71) 

where (a x ) is a slow function of time. Within this ap- 
proximation the Fourier transform of a y in Eq. (|70[) is 
dominated by the Fourier transform of smitA/K) that 
equals 



m[8{u: + A/h) - 6{lu - A/h)} 
so that the integral (|70|) becomes 



/ ; r"77 /A \ 

2^Vl-K) 2 cos^J 



(72) 



(73) 



where T is given by Eq. Substituting this result 

into Eq. (|54[) . taking into account the first of Eqs. (JTTJ) , 
and averaging the resulting equation over fast oscilla- 
tions, (cos 2 (tA/h)) = 1/2, one obtains 



dt 



= r Q (1 - (a x ) 2 ) 



(74) 



This leads to the following relaxation law at t > after 
the system was prepared in the state with arbitrary <r x — 
tanh(roto) < 1 at the moment of time t = 0: 



(a x ) = tanh[r (i + t )} 
sin (ft) 

y cosh[r (t + t )] 

cos (A t ) 

z cosh[r (t + t )] ' 



(75) 
(76) 

(77) 



Our previous consideration of small oscillations of a y>z 
(that is, precession around a x — > 1) corresponds to the 
choice of T to ^> 1, in which case the decay of the oscilla- 
tions is always exponential with the rate To, as has been 
previously found. If the system is prepared in the state 



with a z = 1 (that corresponds to the choice of to = in 
the above equations), it exhibits exponential relaxation, 



2e 



-r t 



cos(tA/h) , 



(78) 



only at Tot 1. The initial relaxation at Tot 1 is 
slower: 



cos(tA/h) 

1w 



1 



(79) 



This later result for a two-state system should be taken 
with a grain of salt, though, as it is likely to be the con- 
sequence of the approximation in which the expectation 
value of the second term in Eq. (|4"5| is replaced by the 
product of expectation values of a y and phonon field. 
Such approximation neglects quantum correlations be- 
tween spin 1/2 and the boson field. In this connection, it 
is interesting to notice that our model can be easily ex- 
tended to a system of more than one flux qubit if all the 
qubits have the same resonance frequency, uj = A/h, and 
are located within a distance from each other that is small 
compared to the wavelength of sound of frequency u>. In- 
deed, for such a system er/2 in Eq. (|4"2"]) gets replaced with 
the total effective spin S = <t 1 /2 + <T2/2 + <T3/2 + .... Since 
the resulting Hamiltonian is linear on S, it commutes 
with S 2 . Consequently, when the number of qubits, N, 
is large, S must behave as a classical large spin of con- 
stant length. In this case, the approximation that ne- 
glects quantum correlations must be good. It leads to the 
same equations (|51[) - (|56[) in which cr is replaced with 
Ner. This amplifies the amplitude of sound by a factor N . 
Consequently, Tq is amplified by a factor N 2 . One im- 
mediately recognizes Dicke superradiancei^ in this effect. 
We, therefore, expect Eqs. §T5§ - (J77]) with T N 2 T 
to correctly describe decoherence in a system of N ^> 1 
closely packed flux qubits. 



Renormalization of the tunnel splitting by the 
elastic environment 



The above consideration shows that decoherence of the 
flux qubit in the elastic environment is dominated by 
phonons of energy A. Meantime, even at T = there 
are zero-point oscillations of the solid that produce elas- 
tic twists. Such twists interact with the flux qubit and, as 
we shall see below, renormalize the tunnel splitting. This 
problem cannot be treated semiclassically as it requires 
consideration of the entanglement of the qubit with the 
excitation modes of the solid. It is based upon compu- 
tation of the quantum average of the Hamiltonian (|42l) , 
(0|/f|0) over the ground state of the solid, |0), that has 
no real phonons. 

Noticing that 



/Q|gJ€kx(oicX-<»L)| 
1, 



kA 



| e -i«kA(akA-Qk A )|0) 

-l ; «-l 2 / 2 , (80) 
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one obtains 



H a = (0|H|0) = - 




(81) 

that is, 

H a = -^a Xl (82) 

where 




is the tunnel splitting renormalized by zero-point quan- 
tum elastic twists. Here £kA is given by Eq. (|4*3")l . 

The sum over k in Eq. (|83p can be computed by re- 
placing it with the integral V J d 3 fc/(27r) 3 . For the two 
transverse phonon modes k x = ±fcekt 2 . Averaging 
over the angles then gives ([ekt]^) = 1/3. Integrating 
over k from zero to fc ma x determined by the size of the 
flux qubit, one obtains 

A eff = Aexp -^>f (84) 

A quick estimate (see Section IIVI) shows that the ex- 
ponent in Eq. (|84[) is always small, thus providing neg- 
ligible renormalization of the tunnel splitting in a flux 
qubit. However, the above result illustrates an impor- 
tant point. If, for some reason, the shear modulus of the 
solid, G = pv 2 , disappeared, this, according to Eq. ([54"]). 
would lead to the disappearence of the tunnel splitting 
as well. The latter is a consequence of the conservation 
of angular momentum: The current cannot reverse direc- 
tion if it cannot transfer momentum to the body. As is 
discussed in the next Section this effect may, in principle, 
be observed in some two-state systems. 

IV. DISCUSSION AND CONCLUSIONS 

We have studied two models that take into account 
mechanical effects associated with quantum oscillations 
of a superconducting current in a flux qubit. These ef- 
fects have simple physical origin. To change direction, 
the current must transfer momentum to the underlying 
crystal lattice. For the current oscillating in a SQUID 
loop, it is a microscopic analogue of the Einstein - de 
Haas effect: The change in the angular momentum of 
the current associated with its magnetic moment must 
be compensated by the change in the angular momen- 
tum of the body containing the current. This inevitably 
entangles quantum states of a flux qubit with quantum 
states of a macroscopic body containing the qubit. One 
can naively imagine that almost instantaneous decoher- 
ence of quantum states of the macroscopic body would 
have a detrimental effect on the decoherence of the flux 



qubit. We show that this is not the case due to the 
selection rule originating from conservation of angular 
momentum. While quantum state of a macroscopic sys- 
tem is, in general, an admixture of a large number of 
rotational states corresponding to different total angular 
momenta, only tunnel splittings of the states belonging 
to the same J contribute to quantum oscillations of the 
superconducting current. Broadening of the tunnel split- 
ting by the rotational states of a qubit is small as long 
as the body is sufficiently large. 

In the first part of the paper we have studied an exactly 
solvable model of a flux qubit entangled with a rigid me- 
chanical rotator. We show that decoherence in such a sys- 
tem is weak due to inverse proportionality of the decoher- 
ence rate, T r = (4hl 2 /T)(kBT/ A), to the moment of in- 
ertia of the rotator, /. To put things in perspective, con- 
sider, e.g., a micron-size flux qubit embedded in a body 
of a comparable small size that is free to rotate. Sound 
of frequency uj — A/H ~ 10 10 s _1 would have a wave- 
length comparable to the size of the body. Consequently, 
in reaction to the oscillations of the superconducting cur- 
rent, such a system would rotate as a whole, making the 
rigid-body approximation developed in Section [XT] a rea- 
sonably good one. Typical value of the moment of inertia 
of a micron-size body is in the ballpark of 10~ 19 g-cm 2 . 
Taking I ~ 10 5 for a micron-size current loop, one ob- 
tains the following values of the parameters in equations 
dS5]) - (J39l): A « 2 x 1CT 8 , p T « 3 x 1(T 4 {k B T / A) 1 / 2 , 
(Too ~ 4 x 10 _8 (fcsT/A). Decoherence is dominated 
by J ~ 10 9 (/cbT/A) 1 / 2 , which corresponds to frequen- 
cies of the rotational Brownian motion ui = hJ/I ~ 
10(fc B T/A) 1 / 2 s" 1 . This provides T ~ 500s" 1 that cor- 
responds to a rather high quality factor of quantum os- 
cillations, Q = A/(hT) - 2 x 10 7 [A/(fc B T)], even in the 
extreme case of a micron size system. 

In the second part of the paper we have studied in- 
teraction of the flux qubit with the twists of the elastic 
body, dictated by the conservation of angular momen- 
tum. Such model has no free parameters. While its 
exact quantum states are not known, one can develop 
a reasonably good approximation in which the internal 
torque produced inside the body by the oscillating cur- 
rent is treated as a source of elastic shear waves. If the 
elastic environment is considered to be infinite in space, 
this is an open system as compared to the closed system 
that consists of a finite-size rotator with a flux qubit. In 
the infinite elastic system the shear waves generated by 
the point source of torque escape to infinity, thus allow- 
ing finite decoherence at T = as compared to the closed 
system. The corresponding decoherence rate is given 
by T e = ; 2 A 5 /(127rfiVt)coth[A/(2fc B T)]. At I ~ 10 5 , 
ui = A/h ~ 10 10 s _1 it is of the order of 10 6 s~\ which 
provides Q = A/(HT) ~ 10 4 . This shows that the effect 
studied in this paper, while allowing weak decoherence, 
can hardly be ignored in designing flux qubits. 

A good check of the validity of the above results can 
be obtained by comparing decoherence rates obtained 
within the rigid-rotator model and within the elastic 
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model. At fc^T > A the ratio of the two rates is 
r e /r r = (47r 4 /3)(//pA 5 ) where A = 2nhv t /A is the wave- 
length of shear waves of frequency oj — A/H. Noticing 
that the moment of inertia of a rigid body of radius R is of 
order pR 5 , we see that T e /T r ~ 1 at A ~ 2R. This agree- 
ment between the two models that consider the same 
effect from two very different angles is quite remarkable. 

In our consideration of the conservation of angular mo- 
mentum, certain effects that may exist in real systems 
have been left out. Among them are interactions of the 
flux qubit with magnetic atoms and nuclear spins that 
can, in principle, absorb some part of the angular mo- 
mentum of the SQUID. For I 3> 1 such processes must 
be suppressed, however, as they require coherent partici- 
pation of many magnetic atoms and many nuclear spins. 
Interaction of the flux qubit with the shear waves of the 
body must be the primary mechanism of the conservation 
of angular momentum. Being unavoidable, it imposes a 
universal upper bound on the quality factor of the qubit. 

The effect of rotations on decoherence can also be un- 
derstood from another angle. At <f> — i*rt that corresponds 
to the uniform rotation of the flux qubit about the Z-axis 
the Hamiltonian (|16[) is equivalent to the Hamiltonian 
of spin 1/2 in the effective magnetic field of amplitude 
A/(2/xb) (a*s being the Bohr magneton) rotating in the 
XY plane at an angular velocity Q = 2lu. Switching to 
the coordinate frame rotating with the field, gives an ef- 
fective constant field applied along the X-axis plus the ef- 
fective bias field in the Z-direction, H'J = —lhioa z — — o x . 
The first term is simply — ui-H\, that appears in the frame 
rotating at the mechanical angular velocity w, projected 
into the | f) an d | I) states. Real bias magnetic field 
B adds the term — B • (/isl) to the Hamiltonian. When 
the field is applied along the Z-axis the full two-state 
Hamiltonian in the rotating (SQUID) frame of reference 
becomes H" = —l(Hw + pbB)o z — ^o x . This proves that 
the rotation of a truncated two-state SQUID system sat- 
isfies Larmor theorem. It is equivalent to the magnetic 
field B/lu = h/fis ~ 10~ 7 Oe/Hz. Effective fields gener- 
ated by slow rotations of the equipment must have neg- 
ligible effect on the flux qubit. However, the effect of lo- 
cal dynamic shear deformations on a microscopic SQUID 
must be noticeable because the corresponding angular 
velocities (ku)(A/H) can easily reach 10 7 Hz, providing 
effective fields in the range of 1G. 

Experiments with flux qubits have shown that signif- 
icant decoherence comes from 1// noise, the origin of 
which has been debate d 15 ' 16 . Notice in this connection 
that relaxation of microscopic shear strains in a solid 
must be a source of dynamical local twists that, according 
to the above discussion, generate local effective magnetic 
fields. It is, therefore, plausible that relaxation of shear 
strains at the location of the qubit is, in fact, responsible 



for the observed 1// noise affecting quantum dynamics 
of the qubit. 

Another observation worth mentioning is amplification 
of decoherence in a system of flux qubits positioned in 
close proximity to each other. This effect may be impor- 
tant in designing architectures of flux qubits if they are 
to be used for quantum computing. It will reveal itself 
when N microscopic qubits with identical tunnel splitting 
A are positioned within the wavelength of sound of fre- 
quency A/H. As has been demonstrated in Section UlIBl 
radiation of sound by such a system and, thus, decoher- 
ence will be amplified by a factor N 2 . This is an acoustic 
analogue of Dicke superradiance that may impose an up- 
per limit on the density of flux qubits. One way to avoid 
this effect in a dense assembly of qubits would be to use 
qubits of significantly different A. 

In Section [III CI we studied renormalization of the tun- 
nel splitting of a flux qubit arising from its interaction 
with zero-point shear deformations. The magnetic mo- 
ment of the current of strength J in a loop of area a is 
ji = Ja/c, which gives I = Ja/(cp,B)- With a = nr 2 
and fc max — 2n/r, the exponent in Eq. (|84l) becomes 
ir 4 HJ I '(3c 2 fi B pv t ) . At J ~ 1/iA it is hopelessly small, 
thus, making this kind of renormalization irrelevant for 
a flux qubit. Notice in this connection that a similar 
effect, described by Eq. (|84|) . may exist for the tunnel 
splitting of the atomic magnetic cluster. In this case I 
would be significantly smaller but fc max would be much 
greater than for a flux qubit. An estimate for, e.g., a 
magnetic molecule frozen in solid He-4 shows that the 
exponent in Eq. (j54")l can easily be of order unity. As 
the He-solid approaches melting transition on decreasing 
pressure, its shear modulus would go to zero, resulting in 
the freezing of tunneling. 

Finally, we would like to notice that treatment devel- 
oped in this paper should apply to nanomechanical de- 
vices incorporating SQUIDs. Such devices have been re- 
cently made and measure d 17 ' 18 . They open the whole 
new field of the entanglement of qubit states with me- 
chanical oscillations. Possible manipulation of supercon- 
ductng qubits by mechanical rotations is another inter- 
esting aspect of the research on nanomechanical super- 
conducting qubits. Our model of a rigid rotator with a 
flux qubit may provide a framework for theoretical stud- 
ies of these effects. 
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